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Abstract 

Introduction: Bone marrow mesenchymal stem cells (BM-MSCs) are multipotent cells that can differentiate into 
different cell lineages and have emerged as a promising tool for cell-targeted therapies and tissue engineering. 
Their use in a therapeutic context requires large-scale in vitro expansion, increasing the probability of genetic and 
epigenetic instabilities. Some evidence shows that an organized program of replicative senescence is triggered in 
human BM-MSCs (hBM-MSCs) on prolonged in vitro expansion that includes alterations in phenotype, 
differentiation potential, telomere length, proliferation rates, global gene-expression patterns, and DNA methylation 
profiles. 

Methods: In this study, we monitored the chromosomal status, the biologic behavior, and the senescence state of 
hBM-MSCs derived from eight healthy donors at different passages during in vitro propagation. For a more 
complete picture, the telomere length was also monitored in five of eight donors, whereas the genomic profile 
was evaluated in three of eight donors by array-comparative genomic hybridization (array-CGH). Finally, an 
epigenomic profile was delineated and compared between early and late passages, by pooling DNA of hBM-MSCs 
from four donors. 

Results: Our data indicate that long-term culture severely affects the characteristics of hBM-MSCs. All the observed 
changes (that is, enlarged morphology, decreased number of cell divisions, random loss of genomic regions, 
telomere shortening) might be regulated by epigenetic modifications. Gene Ontology analysis revealed that 
specific biologic processes of hBM-MSCs are affected by variations in DNA methylation from early to late passages. 

Conclusions: Because we revealed a significant decrease in DNA methylation levels in hBM-MSCs during long-term 
culture, it is very important to unravel how these modifications can influence the biologic features of hBM-MSCs to 
keep track of this organized program and also to clarify the conflicting observations on hBM-MSC malignant 
transformation in the literature. 



Introduction 

Bone marrow mesenchymal stem cells (BM-MSCs) are 
multipotent cells that can differentiate into different cell 
lineages [1]. Human BM-MSCs (hBM-MSCs) are easily 
isolable and are not ethically restricted; thus they have 
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emerged as a promising tool for cell/gene therapy for tis- 
sue regeneration and anticancer treatments. Their applica- 
tion is concurrently tested in various clinical trials [2], but 
their use requires large-scale in vitro expansion, increasing 
the probability of genetic and epigenetic instabilities. 
Spontaneous transformation of mouse BM-MSCs has 
been observed [3-6]; chromosomal instability has also 
been evidenced for rat BM-MSCs [4,7]. Conversely, con- 
founding data exist about the stability of hBM-MSCs and 
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their ability to transform spontaneously in vitro [3,5,8-12]. 
Some authors have reported spontaneous transformation 
of human MSCs, but in several cases, the data were 
retracted, because the occurrence of transformed cells was 
due to cross contamination of the original cell culture 
with tumor cell lines [13-15]. 

Although, to date, hBM-MSCs appear to be less prone 
to malignant transformation during in vitro culture, more- 
detailed studies are urgently needed to evaluate their in 
vitro behavior, particularly as a great variability in terms of 
proliferative capacity and life span was evidenced between 
donors [8]. However, hBM-MSCs have a restricted life 
span in vitro, as does any normal somatic cell, because of 
the phenomenon called the Hayflick limit [16], or replica- 
tive senescence, whereby they exhibit a reduced differen- 
tiation potential, a shortening of the mean telomere 
length, and morphologic alterations [17,18]. It is now evi- 
dent that a strong correlation exists between DNA methy- 
lation-stem cell renewal-differentiation, as well as between 
stem cell culture-copy number changes-spontaneous 
malignant transformation (see reviews [19,20]). Recent 
studies on replicative senescence of hBM-MSCs have 
demonstrated that gene-expression changes are continu- 
ously acquired with increasing passages, influencing their 
differentiation potential [21]. Moreover, DNA methyla- 
tion-pattern variations in hBM-MSCs have been seen to 
overlap in long-term cultures and in aging in vivo, suggest- 
ing that replicative senescence and aging are regulated by 
specific epigenetic modifications [22]. 

The purpose of this study was to track the chromosomal 
status, the biologic behavior, and the senescence state of 
hBM-MSCs derived from eight healthy donors at different 
passages during in vitro propagation. First, we applied the 
conventional cytogenetic technique to observe major 
(>2 Mb) and minor structural abnormalities and to iden- 
tify low mosaic conditions; subsequently, a more-detailed 
whole genomic analysis by array-comparative genomic 
hybridization (a-CGH) was conducted. In addition, the tel- 
omere length was monitored to assess cellular aging in 
vitro. Finally, to evaluate DNA methylation-pattern 
changes after long-term in vitro expansion, a genome- 
wide analysis of DNA methylation was performed compar- 
ing early and late passages, and the results were further 
analyzed by gene ontology (GO) functional analysis. 

Materials and methods 

Isolation, immunophenotyping, and culture of 
hBM-MSCs assay 

hBM-MSCs were obtained from bone marrow in excess 
from eight anonymous healthy donors undergoing mar- 
row harvests for allogenic transplantation at San Gerardo 
Hospital (Monza, Italy). Donor ages were between 20 and 
45 years. An informed written consent was obtained 
from all the subjects, according to the national ethical 



guidelines. Mononuclear cells, obtained after centrifuga- 
tion of the harvested bone marrow in a Ficoll-Hypaque 
column, were suspended in Dulbecco Modified Eagle 
Medium (DMEM; Lonza, Verviers, Belgium) containing 
10% fetal bovine serum defined (FBS; Hyclone, Logan, 
UT, USA), plated in 75 cm 2 (T75) culture flasks, and 
maintained at 37°C in a humidified atmosphere with 5% 
CO2. At this time, cells were considered to be at passage 
0 (PO). After 48 hours, the nonadherent cells were 
removed and the cells attaching to the culture flasks were 
cultured in DMEM plus 10% FBS defined, 2 mAf L-gluta- 
mine, 100 U/ml penicillin, 100 ug/ml streptomycin, 
250 ug/ml fungizone (Lonza), with a change of medium 
every 3 to 4 days. When cultures reached 80% to 90% of 
confluence, cells were washed with Dulbecco Phosphate 
Buffered Saline (PBS; Sigma-Aldrich, St. Louis, MO, 
USA), detached by using 0.25% trypsin in 0.1% EDTA 
(Lonza), and re-iplated (1/3) in 75-cm 2 culture flasks. 

The immunologic characterization of hBM-MSCs was 
performed with flow- cytometric analysis at the Tettamanti 
Foundation laboratories (Monza, Italy) by using specific 
antibodies for the membrane antigens CD33, CD34, 
CD90, CD105, HLA-DR, and HLA-ABC [23]. Osteogenic, 
adipogenic, and chondrogenic differentiation of hBM- 
MSCs was performed by using standard protocols adopted 
in our laboratory [7,24,25]. 

hBM-MSC growth curves 

hBM-MSC growth curves were obtained by plating cells 
on 60-mm-diameter dishes, at a density of 70 to 100 x 10 4 
cells, and counting cells after 24, 48, 72, and 96 hours 
from plating. Population-doubling time (PDT) was calcu- 
lated with the following formula: 

PDT = (CT x ln2)/ln (Nf/Ni) 

where CT is cell- culture time, Ni is the initial number 
of cells, and Nf is the final number of cells. 

Senescence p-galactosidase staining assay 

The amount of senescent cells was determined in hBM- 
MSC cultures, from eight healthy donors and at differ- 
ent passages during in vitro propagation, by using the 
Senescence P-Galactosidase Staining Kit (Cell Signaling 
Technology, Danvers, MA, USA), according to the man- 
ufacturer's instructions. Cells were seeded at a density 
of 1 x 10 4 /cm 2 into 24-well plates and cultured for 24 
hours before senescence-associated P-galactosidase 
staining. The SHSY-5Y neuroblastoma cell line was used 
as negative control of the P-galactosidase staining. At 
the end of the staining procedure, representative images 
were taken from diverse areas of each cell culture by 
using phase-contrast microscopy. For the determination 
of the senescent cell percentage in each culture, an aver- 
age value was calculated by counting, in eight random 
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fields, the total number of cells and the number of cells 
with intracellular blue dye. 

Karyotype analysis of hBM-MSCs 

Metaphase chromosome spreads were prepared from 
70% to 80% confluent cultures at the designated passage, 
according to standard procedures. Actively dividing cells 
were treated with colcemid overnight at 37°C. Cells 
were combined in 1% sodium citrate:0.56% KC1 (1:1) for 
20 minutes at 37°C. Cells were fixed with methanol/ 
acetic acid (3:1). Chromosome analysis was carried out 
by applying Q-bands by fluorescence using quinacrine 
(QFQ banding), according to routine procedures, follow- 
ing the guidelines of the International System for Chro- 
mosome Nomenclature 2009 (ISCN 2009) [26]. On 
average, 25 metaphases were evaluated. 

DNA isolation 

For telomere-length assays, a-CGH, and methylation stu- 
dies, genomic DNA of hBM-MSCs was extracted at differ- 
ent culture passages by using a Wizard Genomic DNA 
Purification Kit (Promega, Mannheim, Germany), accord- 
ing to the manufacturer's instructions. DNA concentration 
was determined on a NanoDrop ND-1000 spectrophot- 
ometer (NanoDrop Technologies, Berlin, Germany). In 
only three cases (P4 of donor 5, P0 of donor 6, and P4 of 
donor 8), the genomic DNA was amplified before array- 
CGH experiment, to enrich its amount, by using a Geno- 
mePlex Complete Whole Genome Amplification (WGA) 
kit (Sigma) according to the manufacturer's protocol. This 
kit allows an approximately 500-fold amplification of 
genomic DNA to be generated. In brief, 50 ng of genomic 
DNA was fragmented and converted to PCR-amplifiable 
OmniPlex Library molecules flanked by universal priming 
sites and then PCR-amplified by using universal primers 
for 14 cycles. 

hBM-MSCs telomere-length assay 

To determine the hBM-MSCs telomere length, Telo- 
TAGGG Telomere Length Assay (Roche Diagnostics, 
Mannheim, Germany) was used at different culture pas- 
sages (P0, P3, P6, P9, and P12), according to the manufac- 
turer's instructions. The positive control DNA supplied 
with the TeloTAGGG Telomere Length Assay is purified 
genomic DNA from immortal cell lines. A total of 1 to 
1.5 ug genomic DNA was digested with an Hinf I/Rsa I 
mixture for 2 hours at 37°C. The sequence specificity of 
these two enzymes ensures that telomeric DNA and subte- 
lomeric DNA is not cut, whereas nontelomeric DNA is 
digested to low-molecular-weight fragments. After DNA 
digestion, the DNA fragments were separated by 0.8% 
agarose gel electrophoresis and transferred, after being 
denaturated and neutralized, to a positively charged nylon 
membrane (Roche Diagnostics, Mannheim, Germany) by 



Southern blotting. The blotted DNA fragments were 
hybridized to a digoxigenin (DIG)-labeled probe specific 
for telomeric repeats at 42°C for 3 hours. The hybridized 
membrane was washed in a high-stringency buffer and 
incubated with a DIG -specific antibody covalently coupled 
to alkaline phosphatase (AP). After the final wash, AP sub- 
strate (CDP-Star), a highly sensitive chemiluminescent 
substrate, was applied and exposed on x-ray film for 10 to 
20 minutes at 25°C. After exposure of the blot to x-ray 
film, the mean size of the different sample smears was 
compared with the molecular-weight marker. 

Array CGH Analysis 

Genomic copy-number analysis was performed by using 
the Agilent Human Genome CGH Microarray 180 K kit 
(Agilent Technologies, Palo Alto, CA, USA) by following 
the manufacturer's recommendations, with a Genomic 
DNA (Female Promega, Mannheim, Germany) as refer- 
ence. The analysis was performed by using Feature Extrac- 
tion vl0.7 and DNA Analytics v6.5 software (Agilent 
Technologies) applying the ADM2 algorithm with a 
threshold of 5, minimum absolute average log2 ratio in 
called intervals of 0.30, and a minimum of three probes. 
Putative chromosome copy-number changes were defined 
by intervals of three or more adjacent probes and were 
considered to be duplicated or deleted when results 
exceeded the ±0.30 range. All nucleotide positions were 
referred to the Human Reference Sequence (GRCh37) 
Assembly Feb. 2009 hgl9 of UCSC. 

MeDIP-Chip 

Methylated DNA immunoprecipitation and chip hybridi- 
zation were performed by following the guidelines of 
Agilent Microarray Analysis of Methylated DNA Immu- 
noprecipitation Protocol (Version 1.0, Agilent Technolo- 
gies). Methylation analysis was performed on a genomic 
equimolar pool of DNA of hBM-MSCs from four differ- 
ent donors: donor 1 and donor 2 at P3; donor 5 and 
donor 6 at P6 (pool of early passages); and donor 1 at P9, 
donor 2 at P10, donor 5 at P12, and donor 6 at P9 (pool 
of late passages). The two pools were used in two inde- 
pendent experiments as reference samples (labeled with 
Cyanine 3); they were hybridized on a Human CpG 
island array (244 K format by Agilent Technologies), in 
competition with the respective methyl-DNA immuno- 
precipitated fractions (labeled with Cyanine 5). The array 
contained 237,220 probes (45- to 60-mer) representative 
of all 27,639 CpG islands in the human genome, at a den- 
sity of about 1 probe per 100 bp. Data analysis was per- 
formed by using Genomic Workbench 6.5 and according 
to the model described by Straussman [27]. In brief, a 
value of Combined Z-score (P value) was assigned by 
Genomic Workbench 6.5 for each CpG island recognized 
by the probes on the array. For each experiment, a 
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bimodal methylation curve was derived: the probe Z- 
scores for each island were averaged to obtain the Island 
Methylation Score (IMS) on the x axis, whereas the num- 
ber of probes was on the Y axis. We then set numeric 
thresholds for determining the methylation status of each 
island. We calculated the distance between the demethy- 
lated (HI) and methylated (H2) peaks and set the upper 
and lower limits for DNA methylation as ±10% of this 
value from the IMS at the lowest point (L), located 
between the two peaks in the bimodal distribution curve. 
Islands with an IMS above the upper threshold were 
assigned a value of +1 (methylated), whereas islands with 
an IMS below the lower threshold were assigned a value 
of -1 (demethylated). Islands with an IMS between the 
two thresholds were considered undetermined (0). 

Gene ontology analysis 

The gene ontology (GO) analysis was performed by using 
GOstat software [28] based on AmiGO (the Gene Ontol- 
ogy database, version 1.8) to identify possible enrichment 
of functional groups, related to "biologic process," in a spe- 
cific input list of genes. The input list may contain genes 
delineated within gain and loss regions detected by a- 
CGH or genes resulted differentially methylated by 
MeDIP-CGI-array. GOstat software output file is a list of 
the P value for each GO term, estimating the probability 
that the observed counts could have occurred by chance. 
GO analysis was selected for the biologic processes, and a 
P < 0.05 was imposed. To limit the number of GO terms, 
a class should comprise more than five genes to be consid- 
ered for further analysis [29]. In addition, GO terms were 
divided into seven functional categories: development and 
differentiation, metabolic process, cell cycle and growth, 
cell signaling, apoptosis and cell death, gene expression, 
and response to stimulus. Categories were ranked in order 
of the percentage of genes found, as described by Liu et al. 
[30]. The percentages of demethylated and methylated 
gene promoters were calculated for each category as fol- 
lows: (total number of genes within a specific category)/ 
(number of genes in the "input list" associated with a spe- 
cific GO term). Because the same gene may belong to dif- 
ferent GO processes, it was counted only once within a 
specific category. Conversely, because the same gene could 
belong to different categories of biologic processes, the 
sum of the percentages could not be 100%. 

Statistical analysis 

Differences in telomere length among passages were ana- 
lyzed by using one-way analysis of variance (ANOVA). For 
each donor and for each passage, a medium spot in the 
range of telomere length (smear) was calculated, and a 
mean value between donors was calculated for each pas- 
sage. Data were expressed as mean ± SEM. Comparisons 
of mean values among the passages were analyzed by 



using a Tukey multiple-comparison test. A 5% probability 
(P < 0.05) was used as the level of significance. 

To assess differences in CpG islands, methylation 
between the two pooled samples of hBM-MSC (donors 
1, 2, 5, and 6), the percentages of each chromosome, 
"early" versus "late" passages, were compared by using 
the Student t test. Differences were considered statisti- 
cally significant with P < 0.01. 

A x test {P < 0.05) was performed to detect signifi- 
cant differences between GO-category percentages of 
demethylated versus methylated gene promoters, mov- 
ing from early to late passages. 

Prospective ethical approval was not sought; however, 
the Ethical Committee of the University of Milano-Bicocca 
analyzed the article retrospectively and recognized an 
overall correct development of the research, including the 
application of the standard informed-consent procedure in 
force at the San Gerardo Hospital. 

Results 

hBM-MSC Characterization 

The hBM-MSCs from eight anonymous healthy donors 
used in this study were characterized according to the 
criteria established by the International Society for Cel- 
lular Therapy [31]: they were plastic-adherent with a 
fibroblast-like morphology (Figure 1), positive for CD90, 
CD105 (Figure 2A), HLA-ABC, negative for CD33, 
CD34 (Figure 2B) and HLA-DR [23], and able to differ- 
entiate toward osteogenic, adipogenic, and chondrogenic 
lineages (Figure 2C through E). hBM-MSCs from several 
donors showed different proliferative capacities, and at the 
same culture passage, the population duplication time 
(PDT) varied greatly from one donor to another. For 
example, at P9, the PDT ranged from 24 hours (donor 5) 
to 97.75 hours (donor 1). Moreover, for hBM-MSCs from 
the same donor, the PDT increased with increasing pas- 
sage (for example, for donor 2 at P6, it was 24 hours, and 
at Pll, it was 68 hours). From P10 onward, for most of 
donors, cultures were characterized by the presence of 
abundant extra- and intracellular debris (Figure 1). The 
achievement of the senescence phase was variable among 
donors (Table 1), as demonstrated by |3-galactosidase 
staining (Figure 3A through E). Negative control repre- 
sented by SHSY-5Y neuroblastoma cell-line culture was 
characterized by |3-galactosidase-negative cells (Figure 3F). 

hBM-MSCs chromosomal profile 

The chromosomal stability of hBM-MSCs was evaluated 
with conventional cytogenetic analysis: at least three dif- 
ferent passages were evaluated, except for donor 3, for 
whom data were available only for P3 (Table 2). All 
hBM-MSCs were generally characterized by a normal 
karyotype with a variable trend to random chromosome 
losses, most likely due to the technical preparation of 
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Figure 1 Human bone marrow stem cell (hBM-MSC) cultures at P3, P6, P9, and P12. At all culture passages examined, hBM-MSCs 
displayed a fibroblast-like morphology. At P9, some extra- and intracellular debris (arrows) appeared and became more evident at P12. Images 
from donor 6 are shown by way of example. Bars, 100 urn. 



chromosomes; conversely, random chromosome gains 
were much rarer. However, in two cases, the presence of 
clonal aneuploidies was evidenced. In one case (donor 4), 
at least 52% of metaphases at P9 presented trisomy of chro- 
mosome 7, also confirmed by fluorescence in situ hybridi- 
zation (FISH) analysis (see Additional file 1, Figure SI); at 
P12, the same karyotype was found in 50% of the meta- 
phases; moreover, in 11% of cells, a loss of one chromo- 
some x occurred, so the total number of chromosomes was 
46. In the second case (donor 5), two equally represented 
subpopulations were evidenced at P4: a normal one, and a 
second with karyotype: 49,XX,+5,+7,+9. However, further 
analysis at P6, P10, and P12 failed to reveal any clonal 
abnormalities, probably because of in vitro negative selec- 
tion of the aneuploid clone. 

Telomere length in hBM-MSCs 

Telomere length was assessed in hBM-MSCs from five 
donors at PO, P3, P6, P9, and P12. For two donors, 



telomere length was examined at PO: in one case (donor 
5), it was comparable to the positive control, whereas in 
the other (donor 8), longer telomeres, not comparable to 
the positive control, were observed (Figure 4A). Regarding 
the next passages analyzed, no differences in telomere 
length were observed by comparing the same passage 
between different donors (Figure 4A). However, differ- 
ences in telomere length were observed between early pas- 
sages (PO, P3) and later ones (P6, P9, P12). As shown in 
Figure 4B, no significant differences in the mean value of 
telomere length were observed between PO and P3, and 
between P9 and PI 2, whereas significant differences were 
observed between P3 and P6 {P < 0.05) and between the 
P0/P3/P6 group and the P9/P12 group (P < 0.001; 
ANOVA test). 

hBM-MSCs genomic profile 

We performed a detailed genomic study for hBM-MSCs 
derived from donors 5, 6, and 8 at P0, P4, and P9 and 10, 
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Figure 2 Human bone marrow stem cells (hBM-MSCs) characterization. hBM-MSC expression of CD105 (A) and CD34 (B) by flow-cytometric 
analysis. hBM-MSC mesengenic differentiation capability (C-E): Alizarin red staining of osteogenic differentiated hMSCs (C); oil red 0 staining of 
adipogenic differentiated hBM-MSCs (D); safranin 0 staining of chondrogenic differentiated hBM-MSCs (E). Bars, 50 pm (D) and 25 pm (C). 



for which it was possible to isolate an adequate number 
of cells. Molecular karyotyping was analyzed by means of 
Agilent Human Genome CGH Microarray 180 K kit (see 
Methods). The comparison of data from three experi- 
ments for each donor allowed us to discriminate between 
large copy-number variations (CNVs), constitutionally 



Table 1 Senescence (3-galactosidase staining assay on 
hBM-MSCs from eight healthy donors 


Donor 


Passage 


% Senescent cells 
(p-galactosidase + ) 


1 


11 


35 


2 


10 


51.25 


3 


-1 


77.5 


4 


10 


50 


5 


13 


17.5 


6 


16 


82.5 


7 


12 


85 


8 


14 


35 



present in the donor's genome, and true chromosomal 
imbalances; Table SI shows this in more detail (see Addi- 
tional file 2). Array-CGH profiles of three experiments at 
different passages from the same donor were almost 
overlapping, and no significant deletions or duplications 
were demonstrated (Figure 5). The overall data showed 
that hBM-MSCs expanded in vitro confirmed a general 
stability of the genomic profile. We performed a GO 
annotation analysis to identify any possible enrichment 
of functional groups of genes within regions with gain 
and loss, but no statistically significant results were 
produced. 

hBM-MSCs epigenomic profile 

To determine the effects of several passages in culture on 
DNA-methylation patterns, we analyzed the methylation 
profile of hBM-MSCs at early and late passages. We 
pooled hBM-MSC genomic DNA from four different 
donors (1, 2, 5, and 6) at early (P3 to P6) and late (P9, 
P10, and P12) passages, to average out any possible 
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Figure 3 Human bone marrow stem cell (hBM-MSC) senescence. p-Galactosidase staining (blue) of hBM-MSCs from donor 2 at P10 (A), 
donor 3 at P4 (B), donor 5 at P13 (C), donor 6 at P16 (D), and donor 8 at P14 (E). SHSY-5Y cell line (F) was used as negative control of the p- 
galactosidase staining. Bars, 50 urn. 



interdonor variation in methylation patterns and outline 
a kind of epigenomic signature specific for the two eras. 
Because we were interested in an overall profile of 
methylation and not in a specific genomic region, we 
applied the MeDIP-CGI-array (Agilent), a microarray 
platform that contains almost all of the CpG islands in 
the genome and that is a reliable solution for distinguish- 
ing highly methylated and unmethylated regions [32]. 
The distribution of CpG island methylation scores in 
both early and late passages shows a bimodal pattern, as 
previously described in a well-validated study [27]. We 
observed that, at early passages, 61.6% of all CpG islands 
were methylated, whereas 38.4% were unmethylated. 
Conversely, these percentages were reversed at late pas- 
sages: 44.7% and 55.3%, respectively (see Figure S4 in 
Additional file 1). To assess differences in CpG islands, 
methylation between the two pooled hBM-MSC samples, 
"early" versus "late" passages, the percentages of each 
chromosome were compared by using Student t test (P < 
0.01). Significant differences were observed in both cases. 
The reversal of methylation percentages between early 
and late passages was observed for all chromosomes, 
except for 18, 21, and X (Figure 6). Moreover, for chro- 
mosomes 4, 8, and 13, a reduction was noted in the dif- 
ference between the two percentages. 

To go beyond the identification of a series of indivi- 
dual genes with a changed methylation profile, we iden- 
tified groups of functionally related genes, based on the 
GO system. Although it is well known that the 



methylation of DNA in 5 promoters suppresses gene 
expression, the role of DNA methylation in gene bodies 
is unclear [33]. For this reason, we decided to limit our 
analysis to only those genes with a change in their 
methylation status in CpG island promoters. In this 
way, the list of "demethylated gene promoters" (that is, 
changing from methylated to unmethylated) contains 
1,284 genes; conversely, the list of methylated gene pro- 
moters (that is, changing from unmethylated to methy- 
lated) contains 518 genes. To simplify the interpretation 
of the large amount of data, we added the same criter- 
ion adopted by Aronica and colleagues [29] (that is, a 
class should comprise more than five genes to be con- 
sidered for further analysis); in addition, we arbitrarily 
categorized the genes altered by extended passages into 
seven functional processes [30]. The Additional file 3 
shows in more detail the two panels of GO terms, one 
specific for demethylated gene promoters and the other 
specific for methylated gene promoters (Additional file 
3: Table S2A and B). The percentages of methylated 
gene promoters belonging to two categories, cell signal- 
ing and apoptosis and cell death, were found signifi- 
cantly different after extended passages (Figure 7) (% 2 
test; P < 0.05). 

Finally, to achieve specific lists of GO terms (that is, 
exclusively present in the "demethylated gene promo- 
ters" or in the "methylated gene promoters" group), we 
eliminated all the redundant GO terms, common to 
these two sets (Table 3). 
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Table 2 Conventional cytogenetic analysis on hBM-MSCs from eight healthy donors at several passages in vitro 



Donor 


Passage 


Norm, 
karyotype 

(%) 


Clonal aneupl. 

(%) 


Clonal 
struct, 
abnorm 

(%) 


Random 
loss (%) 


Random 
gain (%) 


Metaphases 
(number) 


Metaphases/Cells 
(number) 


Mitotic 
index (%) 


1 

(46.XX) 


7 


45.5 


0 


0 


54.5 


0 


22 


24/1 038 


2.3 




9 


66.7 


0 


0 


33 


0 


39 


8/1 000 


0.8 




12 a 
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Discussion 

The chromosomal-genomic profile of hBM-MSCs and the 
potential risk of malignant transformation 

Stem cell-based therapy using hBM-MSCs holds pro- 
mise for treating degenerative diseases, cancer, and 
repair of damaged tissues, for which currently no or 
limited therapeutic options exist. Despite the clinical 
potential of stem cell-based therapy, many risk factors 
were recently described as the "risk profile" by Herberts 
and colleagues [19]. Many identified risks derive from 
the requirement of in vitro expansion and/or differentia- 
tion of hBM-MSCs before administration to a patient. 
Cell-culture conditions may change the characteristics 
of BM-MSCs derived from human and from rat [7,34] 
because of intracellular and extracellular influences. In 



addition, every cell division has a small chance of intro- 
ducing deleterious mutations, and mechanisms aimed at 
correcting these alterations may not function adequately 
during in vitro culture, eventually resulting in a tumori- 
genic phenotype. Several studies on MSCs from different 
sources have highlighted how genomic instability could 
lead to spontaneous immortalization and malignant 
transformation. Spontaneous malignant transformation 
of mouse BM-MSCs after long-term in vitro culture has 
been described [3-6]. Some publications have reported 
spontaneous transformation of hMSCs [11,12]. The 
same authors later retracted their data because of con- 
tamination with immortalized cell lines [13-15], but this 
topic is still open, and further studies are urgently 
needed to ensure the long-term safety of hMSCs. 
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Figure S Example of array-CGH profiles of human bone marrow stem cells (hBM-MSCs) at three different passages of culture. Example 
of chromosome 3 from donor 8 in three overlapping experiments (PO, green line; P4, blue line; P9, red line): the profiles were almost 
overlapping (black arrow, a common copy-number variation (CNV), but P9 showed some exclusive CNVs (dotted rectangle). 



The most important aspect of our study is that, for 
the first time to our knowledge, chromosomal, genomic, 
and epigenomic profiles of hBM-MSCs have been simul- 
taneously evaluated and compared at different passages 
during in vitro propagation. In agreement with the data 
of Bernardo et al. [8], chromosomal stability of hBM- 
MSCs was evidenced for six of eight healthy donors; 
moreover, for two donors (4 and 5) clonal aneuploidies 
were found. In our study, as reported by others [35], the 
abnormal karyotype did not persist on prolonged cultur- 
ing, probably because of clonal selection. Only trisomy 
of chromosome 7 in donor 4 seems to resist long peri- 
ods of culture, also confirmed by FISH analysis (see 
Additional file 1, Figure SI). In addition, by study of 
chromosome heteromorphisms, we can rule out with 
sufficient certainty the possibility of contamination with 
other cell lines (see Additional file 1, Figure S2). The 
aneuploidy of human chromosome 7 might be especially 
interesting because we reported a functional trisomy, 
der(6;6), in our recent study on rat BM-MSCs [7], and 
several syntenic regions exist between these two chro- 
mosomes (see Additional file 1, Figure S3). One might 



speculate that some genes controlling the growth and 
division of cells in these syntenic regions may provide a 
selective advantage, and, indeed, changes in the number 
or structure of chromosome 7 occur frequently in 
human cancers. Even our array-CGH data confirmed a 
general stability of the genomic profile of hBM-MSCs. 
Culture-induced copy-number changes and loss of het- 
erozygosity have been reported for human embryonal 
stem cell lines [35], but the clinical relevance with 
regard to tumorigenic potential of these genomic copy- 
number alterations still remains a matter of debate [36]. 
Moreover, stem cells may be considered potential candi- 
dates for malignant transformation, as some similarities 
exist between their features and those of cancer cells 
[36]. According to the cancer stem cell theory, only a 
small fraction of cells within a tumor (the cancer stem 
cells) are capable of tumor initiation, maintenance, and 
spreading [37]. However, despite the similarities between 
somatic stem cells and cancer stem cells, a direct link 
remains to be shown. 

For many therapeutic applications, MSCs are used in 
an allogenic setting that might facilitate the efficient 



Redaelli ef al. Stem Cell Research & Therapy 2012, 3:47 
http://stemcellres.eom/content/3/6/47 



Page 11 of 1 7 



cn 
cu 
O) 

03 

CO 
cn 
cu 
Q. 



CO 

cu 



c 
o 

eg 

CD 

E 



(/) 
CD 
D) 
CU 
CO 
CO 
CU 
CL 

CD 

-t— ' 

cu 



cz 

o 

on 

0) 

E 



80 



60 



40 



20 



80 



60 



40 



20 



Met 
Unmet 



1 2 3 4 5 6 7 8 9 10 1 1 12 13 14 15 16 17 18 19 20 21 22 X Y 





Met 
Unmet 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 X Y 
Figure 6 CpG methylation profile of human bone marrow stem cells (hBM-MSCs) at early and late passages of culture. Histograms with 
percentages of methylation of each chromosome at early passages (top) and late passages (bottom) of hBM-MSCs. Red, methylated; green, 
unmethylated. 



elimination of transformed cells by the immune system 
[38]. Therefore, the risk of hBM-MSCs would be 
restricted to autologous application in which the 
immune system is less efficient in eliminating the trans- 
formed cells [39]. 

The epigenomic profile of hBM-MSCs and the replicative 
senescence process 

The culture expansion of hBM-MSCs is limited, as it is 
for any other normal somatic cell. After a certain number 
of cell divisions, hBM-MSCs enter a senescent state and 
ultimately stop proliferating. This phenomenon, the 



Hayflick limit [16], is a continuous and organized process 
accompanied by far-reaching alterations in phenotype, 
differentiation potential, and global gene-expression 
patterns [21]. 

Several molecular mechanisms have been implicated in 
this phenomenon [40-42]. The progressive shortening of 
the telomeres has been proposed to be the main trigger 
for replicative senescence; because it functions "as an 
internal clock," with every cell division, the number of 
telomere repeats decreases. Progressive telomere short- 
ening has also been demonstrated for MSCs [17,43]. Our 
data confirm these findings because no differences in 
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Figure 7 Gene ontology (GO) analysis by GOstat software 

Histogram with percentages of gene promoters associated with a 
change in the methylation profile in late passages, classified by 
category of biologic process. Black, percentages of demethylated 
gene promoters; gray, percentages of methylated gene promoters. 



telomere length were observed at early passages (P3 to 
P6) in hBM-MSCs from different donors. In contrast, at 
later passages (P9 to P12), telomeres were markedly 
shorter; furthermore, cell cultures were characterized by 
the presence of abundant extra- and intracellular debris 
and by a decreased proliferative capacity, as shown by the 
reduced PDT and the positivity to |3 -galactosidase stain- 
ing, all clear signs of senescence. However, it is still being 
debated whether telomere shortening is really the initiat- 
ing mechanism or whether it is instead effected by repli- 
cative senescence [44-46]. Telomere shortening may 
induce an antiproliferative signal resulting in cellular 
senescence, which, in turn, acts as a defense against can- 
cer development. 

As an alternative, it has been suggested that molecular 
switches (that is, epigenetic modifications) play a central 
role in regulating cellular aging [47,48]. Moreover, a 
direct link has been described between the maintenance 
of heterochromatic domains, such as those of centro- 
meres and telomeres, chromosome-segregation defects, 
and abnormal telomere elongation [49]. Furthermore, a 
recent study evidenced that DNA methylation patterns 
were maintained throughout both long-term culture and 
aging, but highly significant differences were observed at 
specific CpG sites associated with promoter regions, 
especially in homeobox genes and genes involved in cell 
differentiation [22]. Then, the intimate correlation 
between DNA methylation-stem cell renewal-differentia- 
tion, as well as between stem cell culture-copy number 



changes-spontaneous malignant transformation is now 
evident (see reviews [19,20]). 

In our study, we analyzed the DNA methylation levels of 
hBM-MSCs to delineate a kind of methylation signature 
specific for early and late passages. Because a great varia- 
bility in terms of proliferative capacity and life span was 
evidenced between donors [8], with the aim of eliminating 
interindividual differences, we hybridized a pool of hBM- 
MSC genomic DNA from four different donors on a 
Human CpG island array (Agilent). We revealed a signifi- 
cant decrease in CpG island methylation levels of hBM- 
MSCs during long-term culture, in spite of early and late 
stages being quite close. Furthermore, a reversal of CpG 
island methylated and unmethylated percentages, between 
early and late passages, was observed for all chromosomes, 
except for 18, 21, and X. An explanation for these excep- 
tions could be the lowest content of known protein-coding 
genes for 18 and 21; for chromosome X, it might be unba- 
lanced, compared with other chromosomes, because 50% 
of donors in the pool were male. As it is clear that DNA 
methylation is necessary for controlling stem cell prolifera- 
tion and differentiation [20], and a correlation between 
lineage-specific promoter hypermethylation and lack of 
differentiation capacity toward that lineage was found [50], 
further work is needed to investigate whether the changes 
we found on CpG island methylation levels of hBM-MSCs 
during long-term culture can affect their differentiation 
capabilities. Because it is commonly accepted that DNA 
methylation at CpG islands of gene promoters suppresses 
gene expression, we limited our GO analysis to genes with 
a change in the methylation status (from early to late pas- 
sages) in CpG island promoters. Genes may have an 
antagonistic role in a specific biologic process, because 
some genes must be switched on, and others must be 
switched off, depending on whether they positively or 
negatively regulate this process. For this reason, the per- 
centages of methylation and demethylation gene promo- 
ters could not significantly change for most of the biologic 
processes. 

However, to limit the number of GO terms and to 
simplify the interpretation of the large amount of data, 
we applied the criteria of two published works [29,30]. 
Therefore, we considered GO terms that included more 
than five genes; these were arbitrarily divided into seven 
functional categories and were ranked in order of the 
percentage of genes found (see Additional file 3, Table 
S2A and B). A significant difference between the percen- 
tages of demethylated and methylated gene promoters 
was found only for two categories: "cell signaling" and 
"apoptosis and cell death"; indeed, these two groups 
were only in the "demethylated gene promoters in late 
passages" list (Figure 7 and Table S2A in Additional file 
3). This means that these genes were unmethylated at 
early passages and remained so even at late passages. 
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Table 3 Exclusive gene ontology terms of demethylated gene promoters at late passages of hBM-MSCs 

Category GO terms 

Development and GQ0007275 multicellular organismal development 

differentiation 

GQ0048468 cell differentiation/cell development 
GQ0007399 nervous system development 
GO:0032501 multicellular organismal process 

GO:0009887 multicellular organismal development/system development/organ development/organ morphogenesis 

GO:0030182 generation of neurons/neuron differentiation 

GO:0048699 generation of neurons 

GO:0022008 nervous system development/neurogenesis 

GO:0009790 multicellular organismal development/embryonic development 

GO:0048598 multicellular organismal development/embryonic development/embryonic morphogenesis 
GO:0048666 generation of neurons/neuron differentiation/neuron development 
GO:0009880 embryonic pattern specification 

GO:0048646 anatomic structure morphogenesis/anatomic structure formation 
GO:0000904 cell differentiation/cell development/cellular morphogenesis during differentiation 
GO:0031175 generation of neurons/neuron differentiation/neuron development/neurite development 
GO:0045597 positive regulation of cell differentiation 

GO:0048667 generation of neurons/neuron differentiation/neuron development/neuron morphogenesis during 
differentiation 

GO:0048812 generation of neurons/neuron differentiation/neuron development/neuron morphogenesis during 
differentiation/neurite morphogenesis 

GO:0007389 multicellular organismal development/pattern-specification process 

GO:0007417 nervous system development/central nervous system development 

GO:0045165 cell differentiation/cell-fate commitment 
Metabolic process GO:0016071 mRNA metabolic process 

GO:0009056 catabolic process 

GO:0044248 cellular catabolic process 

GO:0043285 biopolymer catabolic process 

GO:0006401 RNA catabolic process 

GO:0007005 mitochondrion organization and biogenesis 

GO:0009057 macromolecule catabolic process 

GO:0007584 response to nutrient 

GO:0016072 rRNA metabolic process 

GO:0031667 response to nutrient levels 

GO:0009310 amine catabolic process 

GO:0030163 protein catabolic process 

GO:0044270 nitrogen compound catabolic process 

GO:0044265 cellular macromolecule catabolic process 

GO:0006839 mitochondrial transport 

GO:0048878 chemical homeostasis 

GO:0009063 amino acid catabolic process 

GO:0045792 negative regulation of cell size 
Cell cycle and growth GO:0042127 regulation of cell proliferation 

GO:0045786 negative regulation of progression through cell cycle 

GO:0007088 regulation of mitosis 

GO:0000082 G,/S transition of mitotic cell cycle 

GO:0051329 interphase of mitotic cell cycle 
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Table 3 Exclusive gene ontology terms of demethylated gene promoters at late passages of hBM-MSCs (Continued) 

GO:0051325 interphase 

GQ0030308 negative regulation of cell growth 
Cell signaling GO:0007242 intracellular signaling cascade 

GO:0007169 transmembrane receptor protein tyrosine kinase signaling pathway 

GO:0048015 phosphoinositide-mediated signaling 

GO:0019932 second-messenger-mediated signaling 

GO:0045859 regulation of protein kinase activity 

GO:0043085 positive regulation of catalytic activity 

GO:0007167 enzyme-linked receptor protein signaling pathway 

GO:0007219 Notch signaling pathway 

GO:0045860 positive regulation of protein kinase activity 

GO:0033674 positive regulation of kinase activity 

GO:0007243 protein kinase cascade 

GO:0007200 G-protein signaling, coupled to IP3 second messenger (phospholipase C activating) 
GO:0000079 regulation of cyclin-dependent protein kinase activity 
GO:0007165 signal transduction 
GO:0043549 regulation of kinase activity 
Apoptosis and cell death GO:0008219 cell death 
GQ0016265 death 
GO:0006915 apoptosis 

GO:0043067 programmed cell death/regulation of programmed cell death 
GO:0043069 negative regulation of programmed cell death 
GO:0042981 regulation of apoptosis 
GO:0043066 negative regulation of apoptosis 
GO:0008632 apoptosis/apoptotic program 
GO:0006916 antiapoptosis 

GO:0043068 positive regulation of programmed cell death 
GO:0012501 programmed cell death 
GO:0012502 induction of programmed cell death 
Gene expression GO:0010467 gene expression 

GO:0006396 RNA processing 
GO:0006397 mRNA processing 
GO:0008380 RNA splicing 
GO:0006364 rRNA processing 

GO:0040029 regulation of gene expression, epigenetic 

GO:0000377 RNA splicing, via transesterification reactions with bulged adenosine as nucleophile 
GO:0000398 nuclear mRNA splicing, via spliceosome 
GO:0000375 RNA splicing, via transesterification reactions 
Response to stimulus GO:0042221 response to chemical stimulus 
GO:0006950 response to stress 
GO:0009605 response to external stimulus 
GO:0009991 response to extracellular stimulus; 
GO:0006955 immune response 
GO:0010035 response to inorganic substance 
GO:0009628 response to abiotic stimulus 



GO:0042060 wound healing 
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Table 3 Exclusive gene ontology terms of demethylated gene promoters at late passages of hBM-MSCs (Continued) 

Exclusive GO terms of methylated gene promoters at LATE passages of hBM-MSCs 
Category GO term 

Development and GO:0016043 cellular component organization and biogenesis 

differentiation 

GO:0008361 cell morphogenesis/regulation of cell size 
GO:0007276 gamete generation 
GO:0007283 spermatogenesis 
GO:0048232 male gamete generation 
Metabolic process GO:0006512 ubiquitin cycle 

GO:0031323 regulation of cellular metabolic process 
GO:0044255 cellular lipid metabolic process 
GO:0006629 lipid metabolic process 
GO:0006631 fatty acid metabolic process 
GO:0032787 monocarboxylic acid metabolic process 
GO:0043412 biopolymer modification 
GO:0006464 protein modification process 
GO:0043687 posttranslational protein modification 



GO:0008610 lipid biosynthetic process 
Cell cycle and growth GO:0000075 cell-cycle checkpoint 
GO:0016049 cell growth 
GO:0040008 growth/regulation of growth 
GO:0001558 regulation of cell growth 
Response to stimulus GO:0006974 response to DNA-damage stimulus 
GO:0006281 DNA repair 

GO, gene ontology. 

These two categories contain genes that have essential 
functions for the viability and functionality of MSCs 
(that is, the Notch signaling pathway, implicated in mul- 
tiple cell-differentiation processes); thus, they should not 
be turned off. 

After all, to achieve specific lists of GO terms (that is, 
exclusively present in the "demethylated gene promo- 
ters" or in the "methylated gene promoters" group), all 
the redundant GO terms common to these two sets 
were eliminated (Table 3). 

For example, the "metabolic process" class in the 
"methylated gene promoters" group included several 
metabolic processes that could be inactivated with 
increasing passages. Among these, the majority were for 
lipid and fatty acid metabolic process (GO:0044255; 
GO:0006629; GO0006631); interestingly, it was reported 
that adipogenic differentiation potential decreases during 
long-term culture [13-15]. Similarly, the GO:0008361 
(cell morphogenesis-regulation of cell size) was an exclu- 
sive GO term included in the "development and differen- 
tiation" class of the methylated gene promoters group; 
this is quite interesting, because the majority of cells 
acquired a large and flat morphology at late passages. 



Conversely, in the "demethylated gene promoters" 
group were listed genes that were methylated (that is, 
potentially turned off) at early passages and have 
become unmethylated (that is, potentially turned on) at 
late passages. Interestingly, numerous exclusive GO 
terms included in the "development and differentiation" 
class were involved in nervous system development, 
neurogenesis and neuron morphogenesis, and neuron 
differentiation (such as GO:0007399; GO:0030182; 
GO:0048699; and GO:0022008). This finding seems in 
agreement with the hypomethylation of the majority of 
the lineage-specific genes in MSCs reported by others 
[50]. Other potentially turned-on genes belong to the 
"gene expression" class; that is, regulation of gene expres- 
sion and mRNA processing (such as GO:0010467; 
GO:0006396; and GO:0040029). Finally, the GO:0006955 
"immune response" is the only one that has been high- 
lighted as being underrepresented in the list of demethy- 
lated gene promoters. It would be interesting to 
investigate the significance of this observation because the 
immune-modulatory functions of MSCs could also change 
during culture expansion as a result of replicative senes- 
cence [21,38]. 
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Conclusions 

In conclusion, our data indicate that long-term culture can 
affect several biologic features of MSCs. As a consequence, 
in a clinical setting, caution should be exerted before using 
hBM-MSCs for clinical applications. In addition, all the 
observed changes (that is, enlarged morphology, decreased 
number of cell divisions, random loss of genomic regions, 
telomere shortening) seem to belong to a definite program 
that is triggered and finely regulated by epigenetic modifi- 
cations. This developmental process could lead to a reduc- 
tion in the multipotent state of MSCs and might lead to 
tumor formation under specific conditions. It is very 
important to unravel further the epigenetic steps involved 
in this organized program during long-term culture of 
hBM-MSCs; thus appropriate tests should be applied to 
ensure the integrity of the genome and epigenome. 

Additional material 
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